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Abstract 

We discuss scalar brane world cosmological perturbations for a 3-brane world in a 
maximally symmetric 5D bulk. We show that Mukoyama's master equations leads, 
for adiabatic perturbations of a perfect fluid on the brane and for scalar field matter 
on the brane, to a well posed problem despite the non local aspect of the boundary 
condition on the brane. We discuss in relation to the wellposedness the way to 
specify initial data in the bulk. 
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1 Introduction 



Brane world models PU 121 E| lead generically to modifications of the gravitational 
interaction which get reflected into cosmology (see e.g. Aside from the fact 

this can potentially lead to cosmological signatures of extra-dimensions, this has also 
been exploited both for phenomenological constructions and to study more abstract 
questions. On this "abstract" side, one can mention works done in the framework 
of the brane induced gravity model of Dvali Gabadadze and Porrati (DGP model in 
the following) This model can serve as a toy (see e.g. jH]) for exploring some 
long standing problems related to massive gravity, such as the van Dam-Veltman- 
Zakharov (vDVZ) discontinuity |7j. It has also the virtue of leading to interesting 
cosmological consequences for the recent universe jUllE], although the phenomeno- 
logical viability of such a scenario is currently subject to a debate 9j. As far as 
observational consequences of brane world cosmology are concerned, there is cur- 
rently no firm predictions of the possible signatures of these scenario in the CMB 
or Large Scale Structures. This is due in part to the fact that cosmological per- 
turbations on brane worlds are very poorly understood despite the large number of 
works dealing with them and some recent progresses jTOj. Brane world cosmological 
perturbations are also interesting to study in relation with the vDVZ discontinuity 
[TT] . It is some aspects of those perturbations we would like to address in this work. 

We will concentrate on the case of a flat 3-brane of codimension one, and scalar 
perturbations (i.e., scalar from a brane observer point of view). One can first use 
a point of view relying on a 4D projection of 5D Einstein's equations on the brane 
|12j . In this approach the influence of the bulk is encoded into additional degrees 
of freedom with respect to the usual 4D ones, the so-called Weyl fluid [TTH ITU ITB^ 
[HE E] • Knowing the time evolution of those Weyl degrees of freedom is however 
required to know the one of standard cosmological perturbations, but not all the 
Weyl d.o.f. have a local evolution equation on the brane. This means that one 
has in general to solve the full 5D problem in the bulk, and provide some initial 
conditions there to know the evolution of the brane perturbations. In the simplest 
case of a maximally symmetric bulk, one can show that the scalar perturbed Einstein 
equations are solved by a single master variable verifying a second order hyperbolic 
partial differential equation (PDE), the master equation ^H]- The matter equation of 
state, in the simplest case of adiabatic perturbations of a perfect fluid, translates into 
a boundary condition on the brane for the master variable ^Hl This boundary 
condition is "non local" in the sense that it contains derivatives with respect to the 
brane cosmic time of the master variable and of its normal derivative with respect to 
the brane 2 . It is thus not of the usual linear combination of Neumann and Dirichlet 

2 We do not find very well suited the terminology "non-local" to describe the boundary condition 
of interest here, since it only contains a finite number of derivatives, and, as shown in this note, 
the problem can be recast into a standard hyperbolic problem. The "non locality" of the problem 
(as seen by a brane observer) only comes from the fact that the brane is embedded into a bulk 
from where is can receive informations, and not from the form of the boundary condition. The 
latter is indeed "local" when one considers vectors or tensors perturbations ^5]- We will however 
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form. We will show in this work that the same holds true for the case where the only 
matter on the brane consists in a scalar field. The main purpose of this note is to 
show that the differential problem, associated with the master equation formalism, is 
well posed in general, despite the non local aspect of the boundary condition on the 
brane. We will also discuss correlatively the initial conditions to be provided. The 
proof is quite straightforward, and is based on a rewriting of the master equation as 
a PDE of order higher than two in such way the boundary condition on the brane 
takes a standard form. We feel however that a clear presentation of such a proof can 
help a numerical implementation of this problem, but also provides a way, different 
from the usual one, to look at the scalar brane world cosmological perturbations 
as deriving from a standard ("local", that is to say with no derivatives along the 
brane) well-posed hyperbolic problem. This matters in particular when addressing 
issues such as stability of brane worlds [21]. Last, the wellposedness does also play 
a role when using cosmological perturbations of DGP model to study the vDVZ 
discontinuity on cosmological backgrounds [TT] . 

The paper is organized as follows. In the next section, we introduce the differ- 
ential problem associated with brane world scalar cosmological perturbations in the 
Master equation formalism. We then show how this problem can be recast in a well 
posed form. 

2 The differential problem associated with brane 
world scalar cosmological perturbations 

2.1 Background cosmological solutions 

We consider models where our usual 4D space-time is the worldvolume of a 3-brane 
in a 5D bulk space-time. The latter, when not perturbed, will be taken to be a 
slice of a maximally symmetric space-time. In addition we will not consider any 
stabilization mechanism, so that the bulk will be considered empty, and moreover 
we will assume a Z2 symmetric brane. This setup is well suited for e.g. the Randall- 
Sundrum II model [2. [RS model in the following) where the (unperturbed) bulk 
is a slice of AdS§ or for the Dvali-Gabadadze-Porrati model [3], where the bulk 
is a slice of 5D Minkowski space-time. It is also the simplest setting that can be 
considered in order to deal in an exact (and controllable) way with homogeneous and 
inhomogeneous cosmology in brane worlds. The line element of the homogeneous 
cosmological solution reads, in a so-called Gaussian Normal Coordinate system (GN 
in the following), 

ds%\ = —n 2 (t, y)dt 2 + a 2 (t, y)5 i jdx l dx^ + dy 2 , (1) 

where the exact form of the functions a and n, which will not be used explicitly 
here, is model dependent, and has been worked out in the literature [221 II]- in 

stick with this terminology to agree with the literature. 
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this coordinate system 3 , the brane is the y = hypersurface, and the induced 
metric on the brane is FLRW, so that x % can be considered as comoving coordinates 
of standard cosmological observers on the brane. One can in addition choose the 
time parametrization so that n = 1 on the brane, and such that t is the standard 
cosmological time. The function a(t,0) = a^(t) is then the scale factor of the 
FLRW metric on the brane. 

As far as homogeneous cosmology is concerned, the whole space-time can be 
thought of as a slice of a 5D maximally symmetric space-time glued to a copy of 
itself along the brane, the latter being a FLRW space-time. This generalizes (see 
e.g. the well known global picture [21] of a domain wall space-time where 

the metric on the wall (i.e., the brane) is de Sitter |25] . We now turn to discuss 
cosmological perturbations. 



2.2 Scalar cosmological perturbations and the master equa- 
tion 

To deal with cosmological perturbations in the models of interest, one can use f3 
HHHEI] the standard scalar, vector, tensor decomposition (see e.g. |2I3), depending 
on the way the respective perturbations transform under the isometries of the 3D 
maximally symmetric slices of metric (PJ). We will only consider in this article 
the case of scalar cosmological perturbations, which is the most difficult one, but 
also the most interesting as far as phenomenology is concerned. We summarize, in 
this subsection, some results obtained elsewhere on scalar brane world cosmological 
perturbations. 

It has been shown by Mukohyama jTHj (see also ^H]), that the linearized scalar 
Einstein's equations over a maximally symmetric bulk can be conveniently solved 
introducing a master variable Q which obeys a PDE in the bulk, the master equation. 
The latter, when Q has a non-trivial dependence in the comoving coordinates x % and 
when the brane has flat spatial sections, reads in a GN coordinate system (JTJ) 



Vt V M( 5 ) A\ nVL fnQ'Y 



na 3 ) \ 6 a 2 I a 3 \ a 3 



(2) 



where A is defined by A = 5^didj, a dot denotes a derivative with respect to t, a 
prime a derivative with respect to y, and A( 5 ) is the bulk cosmological constant. In 
the rest of this article, we will implicitly consider all the perturbations as Fourier 
transformed with respect to the x*'s, so that the perturbations will all be functions 
of two coordinates (y and t or the characteristic coordinates introduced below, X 

3 Note that this coordinate system is usually breaking down somewhere in the bulk, where the 
functions a or n are vanishing. However the bulk space-time of interest here being non singular, 
it is always possible to choose a better coordinate system for the bulk. Aspects of this issue are 
discussed at the end of the next subsection. 
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and Y) and a comoving momentum k. At the linear level of this work, the different 
modes do not mix, and one can do a mode by mode analysis. In particular, replacing 
A by — k 2 in equation (j2J), the master equation takes the form of a second order (—k 2 
dependent) hyperbolic differential operator acting on y and t dependent functions. 
In the following, we will also use characteristic coordinates X and Y for this operator, 
such that the master equation (J2J) reads 

dxdyfl = U A (X, Y)d x Q + V A (X, Y)d Y Q + W A (X, Y)Q, (3) 

where the coefficients Ua, Va, Wa and the unknown function Q are functions of X 
and Y (and also of the wave number k). For example, when the bulk is a slice of 
5D Minkowski, the master equation reads in characteristic coordinates (X, Y) 



9xW = ^X + AX 2 - (4) 
In those coordinates, the bulk metric reads 

ds 2 = -dXdY + X 2 dx i dx j 5 ij , (5) 

where X is given by a(t, y) and Y is a function of the Gaussian normal coordinates 
y and t which can be easily obtained e.g. by the coordinate change given in |23j . 

The (background) brane trajectory gives a boundary, C(p), in the characteristic 
coordinates (X, Y) plane of the domain, T>, where we want to solve the master 
equation. Indeed, as recalled in the following section, the master variable obeys a 
boundary condition on Cq,) , obtained from the brane matter equation of state. The 
boundary of V is thus defined by C(&) but also by some curve, denoted Cm, where 
we need to specify initial data. This curve can either be nowhere characteristic, but 
also have some characteristic part. We call O the intersection point of Cq,) and Cm 
in the (X,Y) plane (see figure HJ). 

The curve Cq,) is defined by two functions Xq,) (t) and Y^) (t) (t being the cosmo- 
logical time of a standard cosmological observer on the brane) that can be obtained 
e.g. from the relation between characteristic coordinates and Gaussian Normal co- 
ordinate as well as from the cosmological solutions referred to in subsection 12. II The 
exact form of those functions depends on the model considered and in particular on 
the matter content of the brane through the brane Friedmann's equations. However, 
as far as the conclusions of this paper are concerned, the only important point is 
that those functions can be both taken as strictly increasing functions of the cosmic 
time, and in addition, are such that the brane trajectory is nowhere (except possibly 
at the Big Bang, see below) characteristic (meaning here that the brane trajectory is 
nowhere tangent to one of characteristic direction). This holds true for trajectories 
of phenomenological interest. For example, for a brane in Minkowski bulk, as is the 
case for the Dvali-Gabadadze-Porrati model, the brane trajectory is given by 
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Figure 1: Schematic representation of the bulk space-time with the brane trajectory, 
C(6), in the characteristic (X, Y) coordinates. The gray (cyan) part is cutoff and the 
complete bulk space-time is made out of the remaining part, glued to a copy of itself 
along the brane. If a given initial event, 0, is chosen along the brane from which 
one wishes to evolve cosmological perturbations, one needs to specify a boundary 
condition along the brane, and initial data in the bulk along an initial curve C(iy 
We will discuss two different cases, first the case where the initial curve is non char- 
acteristic (left figure), second the case where the initial curve is characteristic (right 
figure). When initial data are specified along [OE[ (resp. [OC[), and a boundary 
condition is supplied along the brane, the solution for the master equation can be 
determined in the (hatched) domain of determinacy T>. 
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It is then obvious to verify that this trajectory obeys the conditions enumerated 
above as long as one has a^) > 0. Moreover, if one assumes that the scale factor on 
the brane, ay,), behaves as t a close to t = 0, one finds that (for < a < 2, which 
is the case for e.g. radiation or matter dominated standard or brane cosmology) Y 
converges as t goes to zero, so that the Big-Bang is indeed in this case a point in the 
two dimensional (X, Y) plane. The coordinate system (jSJ) however breaks down in 
X = 0. A closer look at the way the brane is embedded in the bulk, indeed shows 
|29( I20J |2H] that the brane is a conoidal 4D surface with a summit, an event in 5D 
Minkowski space-time, which is the Big Bang. The brane is also tangent to the 
light cone of the Big Bang along a line which corresponds to the X = coordinate 
singularity. However, one can also go to any point along this line, staying on the 
brane, when t goes to zero, so that the Big Bang also extends along this line. A 
similar picture holds true for the case of an AdS$ bulk (see POJ for a more extended 
discussion of this and of the case where the brane has no Big Bang). 



2.3 Boundary condition on the brane from matter equation 
of state 

So far we have only discussed the bulk perturbed Einstein's equations, as well as the 
background (unperturbed) brane trajectory. We now turn to perturbations of the 
brane induced metric, as well as those of brane localized matter. We discuss how 
to obtain a boundary condition on the brane, from the matter equation of state in 
simple cases of interest. After an appropriate gauge choice, the linearized induced 
metric on the brane can be put in the 4D longitudinal form, which reads (keeping 
only the scalar part) 

ds 2 = -(1 + 2$)rft 2 + a 2 b) (t)(l - 2^)5 ij dx i dx j , (6) 

where dry is the background scale factor on the brane, $ and \l/ are metric pertur- 
bations. The scalar perturbations, ST^ of the matter energy-momentum tensor T£, 
are decomposed as 

5T° = -5p, (7) 
ST? = dM, (8) 
ST] = 5P5) + (s^dkdj - ^5}A\ Stt. (9) 

From the point of view of an observer located on the brane, the functions $, 
Sp, 5q, 8P, 5tc are those which time evolution she would like to know to compare 
with standard 4D results. These functions are all expressible in terms of the master 
variable Q and the anisotropic stress Sir. If one further assumes that the anisotropic 
stress of matter vanishes (this would hold true for matter on the brane consisting 
only in a perfect fluid or a scalar field, those two cases will be discussed in the 
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following), one gets expressions for all the induced metric and matter perturbations 
on the brane in terms of the master variable Q in the form [2U] 



_ 3> _ <E> 

i_5 (o) s_s (i) 



E *(M^ + E Nf M ) W (io) 

with similar expressions for ty,5p,5q,5P. The integers sm and S(o), as well as the 
coefficients K are model dependent. The latter coefficients are functions of the brane 
background induced metric (thus of the brane cosmic time) and the 3-momentum 
square k 2 . They have been given in [201 f° r ^ and DGP models. Expressions (110)1 
are not per se giving a boundary condition on the brane for Q, since their left hand 
side is not known as a function of time. However the sought for boundary condition 
can be obtained from the matter "equation of state" , as we now explain. 

The simplest case one can imagine is the case of an adiabatic equation of state, 
for which one has 

SP = c 2 s 6p, (11) 

where cs is the sound speed. This relation translates into a boundary condition 
on the brane for the master variable O once one replaces 5P and bp in the above 
equation by their expressions as a function of Q and Q' of the form ()10|). This 
boundary condition can be written as 

n (n)+n 1 (n') = o, (12) 

where n and II! are polynomials in cosmic time derivative d t , with cosmic time (and 
A; 2 ) dependent coefficients. Their expression can be deduced from the expressions of 
K p and K p given in [20]. For the RS and DGP models, the degrees of Ho and IT are 
strictly larger than 1, so that the boundary condition (JT2]) is not of the usual linear 
combination of Neumann and Dirichlet form. 

Another example leading to a similar boundary condition is the case where the 
only matter present on the brane is a scalar field. In this case the anisotropic stress 
Stt also vanishes, but equation (Jllj) no longer holds true. Instead one now has 

Sp-6P-5q\2C + 6H\ =0, (13) 

where 5p, 5P and Sq are defined as in equations (J3E1), with being the energy 
momentum tensor of the (canonically normalized) scalar field (with an arbitrary 
potential), and cj) B the background value of the scalar field. This equation can easily 
be obtained from the form of the energy momentum tensor of the scalar field, as 
well as from the conservation equation of its energy momentum tensor. It does not 
require knowledge of the gravitational dynamics (i.e., 4D Einstein's equation need 
not be valid) and can be considered as an equation of state for the scalar field. 
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Following the same path as above and replacing 5p, 8P and 8q by their expressions 
as functions of the master variable Q, one gets a boundary condition on the brane 
of a similar form as equation (|12|) . 

Note that it is only in very special cases, such as the ones studied here, that a 
boundary condition of the form (|T2"j) can be found from the brane matter equations 
of motion (or equation of state). In the most general cases, one do not expect to 
be able to extract from those equations a simple boundary condition for the master 
variable. One should instead solve simultaneously (i.e. numerically) those equations 
and the hyperbolic bulk problem. 

In the next section, we show that despite the non local form of the boundary 
condition (jl2J) . the differential problem associated with brane world scalar cosmo- 
logical perturbations can be recast into a standard hyperbolic form (including the 
boundary condition), and is well posed once one provides suitable initial data in the 
bulk. 

3 Wellposedness 

3.1 Bulk equation and brane boundary condition 

The peculiar features of the boundary condition are that it involves derivatives 
of the master variable of order strictly larger than one as well as derivatives of it 
along the boundary curve, while the PDE in the bulk is of order two. This holds for 
both cases mentioned above, namely boundary conditions derived from equations of 
state (fTTj) and (JT3J). In the RS model, indeed the highest derivative of Q appearing 
in the expression of 5p, 5q and 5P is Q"' [201, which appears in SP. This terms 
appears with a coefficient N(f 2 ) = ( K (5) a (&))~ X (where k? 5 x is the inverse cube of the 
5D reduced Planck mass). This coefficient never vanishes, while all other coefficients 
^0v0> ^(n) an d N(rs) vanish for r + s > 3 (as can be seen from expressions given in 
20 ). Furthermore, one can decompose Cl"' as 

fi"' = (d t Xf (d y X) d 3 x Q + (d t Y) 2 (d y Y) d^Q + ... (14) 

where the dots mean terms proportional to partial derivatives d x dyQ, with r < 2, 
s < 2, and r + s < 3. Because the brane trajectory is nowhere characteristic, d t X 
and d y X do not vanish on the brane (and similarly d t Y and d y Y). One can deduce 
from the above discussion that one can rewrite (|12|). for the RS model, as 

o r r*n= (*(r,s)d r x dpn, (15) 

r+s<r max 

where r max = 3, a(r max ,o) vanishes by definition and the others oc/ r>s ) are known 
functions of the parameter t along the brane (the cosmic time). A similar rewriting 
(with r max = 4) is possible in the case of the DGP model. In this case, the derivative 
of Q of highest order appearing in the expression of 8P is Q" ", and there is no other 
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non vanishing fourth order derivatives in 5p, 5q and 5P. The coefficient N((f 4 ) is 
given by 



2a (6) ft (4) ^1 + — ^r- 



In the above expression k 2 ^, is the inverse square of the 4D reduced Planck mass, 
T is given by T = e/Hr c , where e = ±1 depending on the branch of solutions 
considered for the background, and r c is a distance scale parameter characterizing 
the transition from 4D to 5D gravity in the DGP model (see [201 f° r more details). 
The coefficient N(jf 4 ) does not vanish for physically interesting cases, so that (fT2*|) 
can indeed be recast into the form (fT3j) (with r max = 4). The form of the boundary 
condition (|T5)l leads us to define the (r max + l)(r max + 2)/2 unknowns M( rjS ), with 
< r < r max , < s < r max and r + s < r max by 

«(r,.) = (17) 
and to recast the differential equation Q in the following linear system (for r max > 1) 

dxU( r ,s) = U(r+x j8 ), (18) 

for r + s < r max , 
d x u {r , s ) = dxdy' 1 (U A U( h0 ) + V A u {0 ,i) + W A u {0 ,o)) , (19) 
for r + s = r max and s > 1, and 
9y«(r max ,o) = (^(1,0) + Va«( ,i) + W A u m ) ■ (20) 



Note that to obtain equations (fT9*j) and (j2TJj) . we have derived the master equation 
(jSJ) with respect to the characteristic coordinates X and Y, so that, out of a single 
second order PDE, we are making up a system of PDEs of order higher than two 
(here written in equations (|18H20|) as a linear first order system, in its characteristic 
normal form, see e.g. |2Zj), for which the boundary condition (|15|) takes a standard 
form. The latter reads indeed 

W (r max ,0) = a (r,s)U(r,s), (21) 

with a( rmaxi o) = 0. The characteristics of the first order linear system (fTBH^U]) are 
one family of X = constant lines, and (r max + 1) (r max + 2)/2 — 1 identical families of 
Y = constant lines. Looking at the domain of interest (see figure H}, one sees that, 
in the vicinity of O, there is only one ingoing characteristic in the domain V, the 
X = constant line, for one boundary condition (|2ip. Standard theorems on linear 
differential systems then insure that the differential problem associated with the 
system (|18H20j) and the boundary condition (|21j). is well posed provided one gives 
initial data of Cauchy type on the initial curve Cm, that is to say the values there 
of all the variables U( r , s )- Under this condition, the problem has a unique solution in 
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the domain T> (limited by the brane and the initial curve), and this solution depends 
continuously on any parameter which initial and boundary data may depend on 4 . 
The specification of the initial values for W( r , s ) requires however some care, as we 
now discuss. 

3.2 Bulk initial data 

Given the fact that the master equation in the bulk is of order two, it is indeed 
not possible to choose freely on C(j) all the values of the W(r, a ) without encountering 
inconsistencies. We will consider separately two different cases, first the case where 
the initial curve Cm is not characteristic in O (that is to say that its tangent in 
O is not parallel to the Y axis), and secondly the case where the initial curve is a 
characteristic line (X = constant). We will ask in both cases what freedom we 
have upon the specification of the initial values of the U( r ,s) on Cm. 

If the initial curve Cm is non characteristic in O, it will remain so by continuity 
in some open neighborhood of O, say [OE[ (see figured)). In this case, one can specify 
freely on [OE[ the value of Q and of its derivative normal to C(j) (those are Cauchy 
data on the Cauchy curve C(j) for the second order differential operator defined by 
the master equation (jSJ)). From this specification, the initial values of M(o,o)> u (i,o) 
and M(o,i) are known on [OE[. Because Q obeys a second order hyperbolic PDE, 
the knowledge of Q and the normal derivative of Q along [OE[ is then enough to 
compute the values of all the higher order derivatives, ii( rjS ) with r + s > 2, along 
[OE[. This follows immediately from the fact that [OE[ is not characteristic (and is 
indeed the definition of a non characteristic line). 

We then turn to the case where the initial curve is parallel to the Y axis in some 
open neighborhood [OC[ of O (see figure [TJ. In this case, as is well known, one 
cannot freely specify both Q and its normal derivative on [OC[. Indeed, the PDE 
(J3J) can be integrated along a X = constant segment on which Q is known, to give 
the normal derivative dx^i, on this segment, provided the value of dx& is known 
at one point along [OC[, e.g. in O. This can be generalized by recursion to higher 
order derivatives. Namely if one assumes that one knows along [OC[ the values of 
all the derivatives of order lower or equal to a given m, with m > 1, one first gets 
from those values, and the differentiation of the master equation (J3J), the values of 
all derivatives of the form d x dy^l with s > 1 and r + s = m + 1. One can however 
freely specify the value at O of the derivative d^ +1 Vt. This provide an initial value 
for an ordinary differential equation along [OC[ for the unknown function <9^ +1 f2 
that one can get by differentiating (J3J) m times with respect to X. The latter can 
then be integrated to yield <9^ +1 fi all along [OC[. This freedom to specify the 

4 For this statement to hold, one has to assume some regularity properties on the coefficients of 
the differential equation L, namely Lipschitz continuity |27| . This would however naturally hold 
true as long as the coefficients are continuously differentiable, and that one restrict the domain T> 
to some compact subdomain. Note further that the wellposedness does not imply the absence of 
singularities in the solution 
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derivatives d x Q in O is reminiscent of the wellposedness of the characteristic initial 
problem, for a second order PDE such as (j2J), where one specifies values of Q along 
two intersecting characteristics, X = constant and Y = constant. To summarize, 
in this case, one can specify freely the value of Q all along the initial segment [OC[, 
as well as the values of the derivatives d x Q, for r > 1 at O. If one does so the values 
of the Ut r ,s) w iH De known along the initial curve [OC[. 

Note that in both cases described above, one should make sure that the chosen 
initial data are consistent with the boundary condition (|21|) . Otherwise the solution 
for Q will have discontinuities in derivatives of order r max (see e.g. |27j). Note also 
that once the initial values of ur rtS ) are chosen on Cm, one can compute from the 
expressions (|10|) the initial values on the brane (that is to say in O) of the quantities 
$, Sp, Sq, and SP of interest for an observer on the brane. 

To summarize and conclude, we have shown here that despite the non local 
aspect of the boundary condition on the brane for scalar cosmological perturbations 
of a perfect fluid with adiabatic perturbations or of a scalar field, the differential 
problem associated with Mukohyama's master equation is well posed and can be 
recast in a standard hyperbolic problem, as far as the boundary condition and PDE 
is concerned. Care has however to be taken on the specification of the initial data 
in the bulk, in order to maintain consistency, as explained hereabove. 
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